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the generalized discrete displacement vector 

the generalized discrete velocity vector 

the generalized discrete acceleration vector 

booster (as subscript and superscript) 

payload (as subscript and superscript) 

interface (as subscript) 

non- interface (as subscript) 

with respect to interface (as overbar) 

applied force vector 

reaction vector 

mass matrix 

stiffness matrix 

unit matrix 

modal coordinates 

modal matrix 

constraint modal matrix (interface included) 

damping matrix 

natural frequency matrix 

damping coefficient matrix 

coupling matrices 

time step size 

parameter 

internal load vector 

stiffness kernel 

geometric compatibility matrix 

interface dimension 

non-interface dimension of booster 

non-interface dimension of payload 



INTRODUCTION 


The United States currently utilizes a rather small family of 
launch vehicles (boosters) to support a varied spectrum of satellite 
and spacecraft programs 1 . These launch vehicles have been care- 
fully designed to accommodate a wide range of payload configurations. 

In general* the payload interfaces with the launch vehicle at a lim- 
ited subset of candidate structural "hard" points at the payload 
launch vehicle separation plane. The latest example in the series is 
the Space Shuttle. 

It is important that any candidate payload be designed to with- 
stand the load environment transmitted to the payload from within the 
shielded payload compartment. Such environments commonly originate 
from a static (steady state) vehicle acceleration, a transient or dy- 
namic event such as rocket motor ignition, or an acoustical environment. 
Very often, it is the transient dynamic response behavior of the pay- 
load that constitutes payload design load profiles; hence, it is impor- 
tant that proper attention be given to the payload transient response 
characteristics as influencing major design decisions. 

Present analytical techniques by which such design loads are pre- 
dicted are very costly and time consuming. The calendar turnaround 
time of a given cycle usually is lengthened when the payload design 
organization, the booster organization and the payload integration 
organization are different companies. Indeed, a fair amount of coor- 
dination is necessary to make the transfer of information between 
those three organizations optimal. Unfortunately, this coordination 
is very difficult to establish, resulting in considerable time delays. 
Moreover, these costs and delays repeat themselves for every load cy- 
cle (i.e., every time a change is made in the booster or payload). A 
typical example is the development of the Viking Orbiter System [2]. 

The ever increasing number of modal coordinates necessary to 
model today's aerospace structures not only increases the cost of 
a load cycle, but also imposes greater demands on the analyst to keep 



the models within range of current computer capabilities. 

The objective of the present work is to develop a "full-scale" 
payload integration method which reduces the cost of a load cycle 
dramatically and will be capable of handling very large systems. 

This new approach is a "full-scale" method in the sense that it ac- 
tually solves the coupled booster/payload system equations and does 
not involve any additional approximations or assumptions as compared 
to the standard transient analysis. 

The present report is part of a research effort into possible 
developments of so-called "short-cut" methods. A "short-cut" me- 
thod essentially introduces certain assumptions and/or approximations 
which partially or completely circumvent the coupling between the 
booster and the payload. As will be shown in Chapter II, it is indeed 
possible to derive such "short-cut" methods. However, a "short-cut" 
method inherently produces approximate system responses and, therefore, 
approximate design loads. The question then is if the increased cost- 
effectiveness of such a short-cut method justifies the loss of accur- 
acy in the response and the design loads. 

The present "full-scale" method is very cost effective and 
directly applicable for the shuttle payload design case. This develop- 
ment effective analysis tool without the use of approximation methods. 
For a comparison of "short-cut" and "full-scale" methods, the reader is 
referred to References [3] and [4] . 
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CHAPTER I: A DIRECT INTEGRATION METHOD FOR LOW FREQUENCY 
ENVIRONMENTS - METHODOLOGY 


1. INTRODUCTION 


The shear size and complexity of today's aerospace 
structures has created the need for better and more cost- 
effective payload integration techniques. Indeed, the 
cost associated with a typical design load cycle, has 
become quite substantial. In order to arrive at the 
optimum design, several of these design load cycles are 
necessary for a particular event. In addition, for each 
event various load conditions must be investigated which 
also add to the total cost of the payload design effort. 

This section derives methodology of a payload integration 
method which significantly reduces analysis and design 
turn around time while retaining the accuracy of a "full- 
scale" method. 

The configuration of multiple payloads connected to the 
booster through separate interfaces is typical for most 
Shuttle missions. The fact that many of these payloads 
are not directly coupled allows for a significant sim- 
plification of the booster/payloads systan equations. 
Superfluous interface degrees of freedom on the booster 
side can be accomodated within the formulation. Super- 
fluous interface degrees of freedom are those interface 
degrees of freedom which are included in the launch vehicle 
model but are not connected to any payload. The superfluous 
interface degrees of freedom arise because the booster 
organization cannot afford to reconstruct a booster model 
every time the interface with payload(s) changes. 

A numerical integration scheme used to obtain the booster/ 
payload(s) system response is defined. The standard ap- 
proach is to obtain the so-called "modal modes" i.e. the 
coupled system modes in order to decouple the system 
equations. The present approach avoids the solution of 
such a system eigenvalue problem. A Newmark-Chan-Beta 
numerical integration scheme is used to directly determine 
the system response. This technique takes advantage of 
the peculiar structure of the equations of motion for the 
system. All zero entries of the system mass and stiffness 
matrices can be eliminated. In fact, a comparison feature 
can be iiplemented so that elenents close to zero can be 
omitteu, reducing the cost of the integration routine even 
more. The comparison feature makes the full-scale method a 
so-called short-cut method. A full-scale method can be 
defined as a method which yield "exact" results compared to 
the standard transient analysis technique, whereas a short- 
cut method introduces an approximation or assumption which 
leads to a more cost-effective but less accurate solution. 

A Fortran computer program has been written and implemented 
on the CDC Cyber 172. The final remarks section discuss 
advantages and disadvantages of the proposed approach. 
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FORMULATION OF THE PROBLEM 


The objective of this section is to describe the payload inte- 
gration problem as it occurs in many of today's engineering applica- 
tions. 


As an example, let us consider the landing of the Shuttle Or - 
biter (» the delta-wing -airplane-like module) carrying a certain 
payload (e.g. the space telescope). Obviously, when the orbiter 
touches the landing strip it will experience reaction force's. These 
forces will be transmitted to the payload through the interface (i.e. 
through the connection points between the orbiter and its payload) . 

The payload then, will undergo elastic displacements. The question 
then is, will the payload be eU3le to withstand those displacements 
without being damaged? The answer to this question requires a dynam- 
ic analysis of the booster/payload system. 

First, an analytical booster and payload model is developed. 

This involves the construction of mass and stiffness matrices, using 
any suitable method (e.g. the finite element method). These models 
are usually extremely large ( i.e. a very large number of degrees of 
freedom) and must be reduced to a suitable working size. The re- 
duction of the booster and payload models requires a considerable 
€unount of engineering judgement. Indeed, the objective is to retain 
enough information or fidelity in the original model and to still arrive 
at an acceptable working size. 

At this point it is possible to derive a coupled booster/payloac 
system model. The key to obtain the system equations is the fact 
that the displacements at the booster side of the interface must be 
the seune as the displacements at the payload side of the interface. 
Similarly, the reactions at the interface must be equal in magnitude 
but opposite in sign. Although the booster/payload system equations 
could now be solved, in most applications this system would still 
be large and costly to solve. Fortunately, in many cases the force 
environment has a rather low frequency content. In our example of the 
orbi* ar landing, this would mean that the reaction forces on the or- 
biter from the landing gear would have a low frequency content (e.g. 



< 50 Hz). It is a K.iown fact that if the forcing function has a 
frequency content below say 50 Hz, that the system response above 
50 Hz will be negligible compared to the response below 50 Hz. This 
at least opens the possibility of truncating the modes and frequencies 
in the booster and payload model in order to further reduce the size 
of the system. Although certain precautions must be taken, this ap- 
proach is indeed a valid one and leads to very accurate results. 
Consequently, the standard procedure is to solve an eigenvalue prob- 
lem for both the booster and the payload and write the system equa- 
tions in terms of these modes and frequencies. The size of the sys- 
tem can then be reduced by truncating all the modes above the so-called 
cut-off frequency (* 50 Hz in our example). As mentioned above, 
certain precautions must be taken in order not to lose essential in- 
forriiation, especially at the interface. 

The next step is to solve the set of system differential equa- 
tions. This step is usually a very expensive one. The standard ap- 
proach is to solve the system eigenvalue problem in order to decouple 
the system. The solution of the uncoupled system equations is then a 
relatively cost-effective process. Although the system modes can 
again be truncated according to the cut-off frequency, it is still 
a very expensive it«n to obtain these sy tern modes. 

Once the system response is calculated, one is now in a position 
to calculate payload member loads. Indeed, Hooke's law defines the 
relationship between the displacements in the member and the loads 
in the member. Note that in practive, one often uses the so-called 
"acceleration" method to obtain loads, as we shall see later in this 
report. 

Finally, the member loads can be used to obtain maximum stresses 
and strains and enable the payload designer to determine whether or 
not the members will be damaged during the landing of the orbiter. 

Of course, several other events must be evaluated (e.g. ignition, sep- 
aration,meco, etc.). Moreover, the response/load problem is an itera- 
tive one, i.e. that once the designer changes the payload based on 
a set of design loads we have a new booster/payload system and, there- 
fore, a new set of laods. In general, several so-called load cycles 
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are necessary to arrive at the optimum design. The flow chart on page 24 
graphically despicts the essential elements in a load cycle. 

From the above, we can conclude that the dyneunic analysis part 
in a design effort is often a very costly itera. Consequently, 
much research has been done in order to improve and make the lo^ cy- 
cle more cost-effective [4^ The present report is such an effort 
The objective is to make the response-link and the loads-link more cost- 
effective without sacrificing any accuracy. 
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DERIVATION OF THE BASIC EQUATIONS 
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The purpose of this section is to derive the basic equations 
needed for a booster/payload integration effort. 

Figure 1 shows the free body diagrams of the booster B and the 
payload P. The booster and the payload are connected to each other 
through the interface. Physically, the interface is the collection 
of structural "hard points" which the booster and the payload have in 
common. Mathematically, this means that 





for all times t 


( 1 ) 


the 


From the free body diagrams in Figure 1, we can easily write 
equations of motion for the booster B and the payload P as. 
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Note that the equations of motion for B and P are still uncoupled 
and still contain the unknown reactions {Rj} and | Rj[ . Several well- 
known procedures [4] exist to derive the coupled set of system equa- 
tions, i.e. to eliminate the unknown reactions and redundant inter- 
face displacements | x ^ | . One such technique uses interface restrain- 
ed booster and payload displacements and as will be pointed out later, 
is ideally suited for our purposes [5]. 


The basic idea is to construct the following transformations. 



( 3 ) 
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OF POOR QUALITY 



( 4 ) 


where 

Cs ] - - Ck^ ] ^ [kJj] (6) 

Using these transformations together with Equation (1 ) , it is pos- 
sible to write 



This final transformation (7) substituted into Equation (2) will 
eliminate the redundant set of displacements | | and in the pro- 

cess it will also eliminate the unknown reactions |Rj| abdjR^ yield 
ing. 
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( 8 ) 


where we assumed that no forces afe acting on the payload P (an as- 
sumption which need not be made) . Equation (8) represents the coupled 
booster/payload system of equations of motion in discrete coordinates, 
using non-interface coordinates for both the payload and the booster. 

Note that damping can always be included. Furthermore, the matrices 
T T 

fTgKgTg] and PpKpTp] are zero only when the interface is determinate. 

B P 

Next, we introduce the interface-restrained modes and 

and also the interface modes [4>j] / i.e. 



1 I ' ! ^ I * I 4 1 
1 -I I - I 1 


( 9 ) 


Using Equation (9) in Equation (8) leads to 
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Note that we also introduced damping. We shall ma^e some more remarks 
about damping in Chapter II. 

At this point one usually truncates the booster and payload modes 
according to a predetermined "cut-off frequency." This cut-off fraquen- 
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cy is based on the significant frequency content of the forcing func- 
tion . Assuming there are n frequencies below the cut-off 

frequency, one usually retainsvf x n modes. This is a rule of thumb 
developed from observation of practical modal coupling cases. This 
rule of thumb allows for a safe margin. Indeed, often times modes 
with frequencies higher than the cut-off frequency couple into 
modes with frequencies smaller than the cut-off frequency and bring 
about significant perturbations. 

Finally, as will be explained in Section 5 of this chapter, 
it is imperative not to truncate interface modes, even though they 
might contain modes with frequencies higher than the cut-off fre- 
quency. 
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Figmre 1 Frte-Body Diagrams of Booster B and Payload P 
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THE nl:'.erical integration scheme 
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The objective of this section is to choose a suitable integra- 
tj .m scheme to directly integrate Equation (10]L This is in contrast 
. th the usual approach, where a new eigenvalue problem with [m] 
t'ld [k] as mass and stiffness matrix is solved. 

It is important to recognize that system (lO)can have a frequen- 
cy content much higher than the cut-off frequency of the forcing 
tt rm. Usually, those higher frequencies will not produce signifi- 
c nt responses. Therefore, a suitable niimerical technique should 
b-.^ capaUale of using a stepsize h, which reflects only the highest 
frequency of interest but at the same time remains numerically sta- 
ble. For example, a Runge-Kutta routine would not be suitable be- 
cause it requires a time-step consistent with the highest frequen- 
cy in the rystem, even though these high frequencies may not be of 
interest to the analyst. Although there are techniques to obtain a 
gocd estimate of the highest system frequency, using such a fre- 
quency to determine the stepsize would unnecessarily increase the 
cost of the response routine. 

A mew..od that satisfies cJaove constraints is given by the New- 
nark -Chan- Beta integration method: 

•I'Sl i ^ 6hV>i*i (15) 

it! *Cc]{i3| +[K])q| * 1*' t i+i ( 15 ) 

2 

where the method is vuiconditionaliy stable if 6 ^ (2a+l) /16. Arti- 
ficial positive damping is introduced when a > 0.5, and artificial 
negafve damping if o < 0.5. Good values for our purposes are o = 0. 
and 13 = 0.25. Theoretically, the time step h can then be given any 
value while the scheme remains stable. In fact for very large values 
of h, the scheme generates the static solution of Equation ( i o ). Also 
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the scheme will damp out the highest (and least in^rtant) modes 
while preserving the lower ones. In addition, as we shall see 
shortly, the Ne%rmark-Chan-Beta scheme is capable of taking advan~ 
tage of the peculiar structure of the present equations of motion 
Indeed, let us substitute Equations (14) and ( 15 ) into (16) and 
obtain 

[M + ohC + 8h^] I I i+1 ■ [C + hK] I <5 I . 

- [(l-O)hC + (J 5 - 8 )h^K] I *q*[ ^ - [K]{q[ ^ 


or, using Equation (11) , we can write 
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are diagonal matrices. Consequently, the evaluation of ^li' 

and are very cost-effective. Also, note that (i=l,...,9) 

are one-time calculations. 

Normally, the solution for | b*| in Equation (13) requires a 

triangular decomposition and must be repeated for every h. However, 
in this case the unique form of the coefficient matrix of 
makes it possible to avoid such decompositions. Moreover, it is 
possible to completely take advantage of the diagonal and zero par- 
titions appearing in that coefficient matrix. 

First, let us premultiply Equation (18) by [-BD”^ I -BD”^] 
which yields the following expression forjqj^^^}. 
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I “ {^1 ^ 2 ^Bi ^li ^3^Bij 

(23) 

with 

^1 = 

[D^ +AjB^ + 

( 24 ) 

^2 = 

-BD^^ , A 3 = -PD 3 ^ 

( 25 ) 

Furthermore , from 

Equation (18), we easily obtain 




(26) 

!*'4i+i! 


(27) 


Equations (23) , (26) , and (27) represent the final set of 
recurrence relations replacing Equation (1B\ Note that Equation (24) 
represents the inversion of an IFxIF matrix where IF is relatively 
small in many applications. This effectively removes the problem 
of triangular decomposition of a large matrix. Also, note that the 
cost of the algorithm primarily comes from multiplications involving 
matrices and A^, Note however, that their dimensions are 

IFxIF , IFxNB and IFxNP respectively. Again, in many cases IF is 
a rather small number. In addition, the routine requires much less 
core memory which allows for the solution of much larger proolems. 
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The purpose of this section is to briefly review the "acceleration* 
approach to calculating loads and at the same time point out some pos- 
sible savings. An elementary member load transformation can be written 
as 

I L I « [k'i'] I Xp( (28) 

Therefore, once the system response is known one can substitute the dis- 
placement vector |Xp| into Equation C2 8) and obtain the member loads. This 
direct approach is called the "displacement" method. In many cases, 
this is a perfectly valid approach especially if all the modes in 
Equations (9) and 0-2) are kept. However, if modes are truncated accor- 
ding to a cut-off frequency this procedure often leads to inaccurate 
results. This can be corrected by using the so-called "acceleration" 
method whereby Xp is replaced in terms of applied forces and acceler- 
ations using the system equations of motion v ‘thout damping. 


Hruda and Jones [8] introduced a load transformation consistent 
with modal synthesis techniques. In terms of the present notation 
Equation (28) can be replaced by 


where 


{Lf 


[LT]3 



[LT2] {x®j 


(29) 



[LTl] = 

[k-i-lClpEplJlC-Mpldpl^i Tp) 

(30) 


[LT2] = 

[kT][Tp] 

(31) 

and [Ep] 

- [IpKpIp 

] Note that [LT2] = 0 

when the interface is 


determinate. 


Normally, the displacement vector in Equation (29) must be 

written in terms of forces and accelerations using the second parti- 
tion in Equation ( 8) . Thi s not only increases the cost of the proce- 
dure, but also introduces the booster accelerations | q^| and forces | Fg 
into the problem. We shall now show that this is not necessary if all 
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interface modes are retained in 



From ’'quation (8 ) , for an indeterminate interface we can write 

!*?1 • ' tV ' {'' b ! 

-[T^Tb + T^MjTj,] jl^j ) 

At this point one can introduce the modal accelerations using 
the transformations { 9) , yielding 

-[^''b'b^n’ I’n ! - (33) 

-[t^m^Tb + TX^p]e«?iht| ) 


This is the expression that ordinarily must be used in Equation 
(29). However, if and only if we keep an the interface modes in ] 
we can write that [TgKgT^ + ^pK^Tp] ^ 

In addition we can then write |x® | = [4>®]|q® | without any loss 

of accuracy. Using these facts in Equation (33 ) and taking damping 
into account, we obtain: 
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which 
( 10 ) . 


is exactly what we ohl 
Therefore, Equations 



e^.n from the second partition 
(29-31) can be replaced by: 


[LTl] 




[LT2]{ q^} 


of Equation 


(35) 
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[LTl] = 

[k'?][lj,Ej,I^][-Mj,](Ij,l^; Tj.4®) 

(36) 

[LT2] = 

[kf][Tp4®] 

(37) 


Next, let us write Equation (3g) as follows: 

[LTl] = [k'i'lCiplc-CEplCijMpiplCi^] j 

T B 

- [Ep][ljMpTp][$®]) 

It is now relatively easy to show that [Ep][lpMpIp][$^] = Ct^][wp]”^ 

even when not all the payload modes are retained in Therefore, 

Equation (38) can be written as. 


[LTl] - 

- [Ej,][^MjTj,][t»]! 


(39) 


Using Equations (35), (37) and (39) has several advantages over 
the approach outlined in Reference [8]* First, in case the interface 
is indeterminate (i.e. if LT2 ^0) it is not necessary to use Equa- 
tion (32 ) for I X j ^ , if and only if we keep all the interface modes 
[$j]. This not only makes the evaluation of Equation (35) much simpler 
but also reduces the amount of information to be stored in the course 
of the response calculations ( i.e. only | q^ | and ^q® • must 

be stored). Secondly, the rather expensive evaluation of [a] = 

-[k«][lp][E 
by tne more 

ing unit loads successively to each of the non-interface dofs. . often 
becomes a rather expensive item, considering the potentially very large 
payload models. The other term in Equation (38 ) » [y] = “[k't'][lp][Ep] 
[lpMpTp][$jJ can easily be evaluated by first forming the product 
[l^MpTp][t®] , the columns of which can be looked upon as inertial loads 
applied at the interface of the payload. The corresponding def].ec- 
tions are equal to [e] = [Ep][lpMpTp][4>j] . Note that [Ep] need not 
to be recomputed, i.e. no decomposition of C^p^p^p] is necessary since 


p][lpMpIp][l^] as proposed in Reference 8 is now replaced 

r . T _ rTPir 2 -,-] Indeed, apply- 


efficient computation of [ 8 ] “ C'l’^][wp] 
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this has already been done in the modes calculations. Finally, one 
can also evaluate [e] = [Tp$“] inexpensively in Equation ( 37 ) so that 
we obtain 


[LT]] = [k'F]([3].[6]) 

(40) 

0L.T2] = [k'l'][e] 

(41) 


Observe that the payload organization can easily save [ 6 ], [5] 
and [e] so that any member in the payload can now be investigated with- 
out recalculating these matrices . Observe that Ref erer.'je 8 requires the 
payload analyst to make the choice of members before the evaluation of 
Equations (40) and (41). If for some reason an additional member has 
to be investigated a reevaluation of LTl and LT2 is necessary. 

Finally, it should also be noted that the loads calculation does 
not involve the "modal modes" which reduces the computational cost 
even more. 


’.iS 


Qf POOR 


- 20 - 



6 . 


CONCLUSIONS AND SUMMARY 


A new payload integration approach ha? been presented. It is 
a "full-scale" approach in the sense that it does not introduce new 
?ssumptions or approximations compared to the conventional "exact" 
solution techniques. Improvements over the conventional techniques 
are introduced in both the response and loads calculations. 

The response analysis uses an adaptation of the Newmark-Chan-Beta 
numerical integration technique. This integration =^rneme is directly 
applied to the coupled system equations (i.e. booster/r oad system) 

thereby avoiding the expensive solution of a system e: Mue prob' 

lem. The Newmark-Chan-Beta technique has the convenien ature that 
the step size can be based on the "cut-off frequency" associated with 
the forcing function regardless of the highest system frequency. This 
particular feature is necessary in the present method because the high- 
est system frequency is not known a priori. Although there are tech- 
niques to determine the highest frequency, it is very likely that this 
highest frequency will be much larger than the cut-off frequency which 
would lead to a much smaller time step. It should also be noted that 
the present approach allows for the solution of much larger systems. 

Next, we derived a load transformation consistent with the above 
modal synthesis method. Several cost saving features were introduced. 
First, we showed that in the case of an indeterminate interface it 
is not necessary to write the interface displacements in terms of accel 
erations and forces, provided one keeps all the interface modes. In- 
cidently, one could actually keep the discrete interface displacements 
instead of introducing the interface modes. Secondly, it was shown 
that several simplifications can be affected in the calculation of [LTl] 
and fiT2] leading to a more efficient and convenient loads calculation. 
Sihally, it should also be pointed out that we do not involve system 
modes which reduces the cost and simplifies the analysis. 

Tt is estimated that the present approach will reduce the computer 
cost of a payload integration effort by a considerable amount. Con- 
sidering the numerous load cases that must be considered in the course 
of a design effort, the present approach may prove to be of great value 
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It should also be pointed out that this technique can easily be 
adapted to a "short-cut" method. Indeed, Section 4 of Chapter II 
will be devoted to outlining several pocsible short-cut procedures. 

Finally, also note that the method as presented in Section 4 is 
considerably more cost-effective than the one presented in Reference 
10 . 
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CHAPTER H i A DIRECT INTEGRATION METHOD FOR LOW FREQUENCY ENVIROH> 
MENTS - IMPLEHENTATION 

1. INTRODUCTION 

This chspttr discussts In gtntral ttms tho softwart packagt 
associatad with a coinplatt boostar/payload rtsponsa and loads 
analv*-1s. An attampt will ba mada to claarly link tha thaory of 
Chaptaf I with tha spacific program and subroutina dascriptlons. 

This will gfva us the opportunity to touch upon soma of tha const- 
raints and difficultlas Invariably associatad with tha davalopmant 
of a practical payload Intagratlon sofbvara packaga. Soma factors 
tc ^onsldar ara: computar ccra usagai convarganca; avallabla data; 
tha saparatlon of boostar, payload and Intagratlon organizations; 
work schadulas; anginaaring tima; aasa of program usaga; computar 
cost and ralctad afficlancy of algorithms; rausa of axisting Infor- 
mation; raquirad accuracy vaAsus cost; handling of potantlally larga 
modils; ate. 

faction 2 of this chap tar prasants a ganaral dascriptlon 
of tha organization and componants of tha softwara packaga. In 
particular, wa axplain tha purposa and contants of tha componants 
and how thay ralata to aach othar. 

Saettun 3 peasants a simpla sampla problam and how It Is 
analyzad and avaluatad. Also, prallmlnary rasults of mora raallstic 
analysas will bo Included. 

Finally, conclusions and possible short-cut approaches ara 
Introduced In Section 4. 
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2. ORGANIZATION - GENERAL DESCRIPTION 

This section outlines the organization of the software package. 

The figure 2 represents a flow diagram of a complete booster/payload 
integration problem. Each of the flow diagram blocks has a program 
associated with it. Therefore, there are six programs: PROGRAM 
BOOSTER, PROGRAM PAYLOAD, PROGRAM INTFACE, PROGRAM FORCE, PROGRAM 
SYSRESP, AND PROGRAM LOADS. Each of these programs draws on a pool 
of subroutines called FORMA (FORTran Matrix Analysis). FORMA is 
a library of subroutine coded in FORTRAN IV for the efficient solution 
of structural dynamics problems. These subroutines are in the form 
of building blocks that can be put together to solve a large variety 
of structural dynamics problems. The FORMA library was developed by 
the Dynamics and Loads Section at Martin Marietta Aerospace and is 
being updated and expanded whenever the need occurs [11] . 

It should be pointed out that other libraries can be used and that 
the proposed integration method does in no way inherently depend on 
the FORMA library. However, in this report, the software is built 
around the FORMA subroutine library (in particular, the Partition- 
Logic version) and therefore, the user is assumed to have a working 
knowledge of that library. • 

There are several reasons motivating the PROGRAM apporach. Be- 
cause all FORMA routines ere written in terms of variable dimensions, 
it is possible to write each PROGRAM for the specific dimensions of 
the problem at hand, thereby optim'' 2 ing computer core usage. Also, the 
user often has at his/her disposal data already generated by other means. 
The PROGRAM approach allows the user to omit the recalculation of 
that particular data. For example, a set of so-called expanded modes 
could be available. The user then can directly read those expanded 
modes into the PROGRAM and omit the use of a subroutine which calculates 
those expanded modes. Furthermore, the PROGRAM approach allows for 
the separation of data generated by the booster, payload and integra- 
tion organizations. Indeed, very often there three organizations are 
physically at different locations and data peculiar to one organization 
often is not readily available to the other organizations, PROGRAM 
BOOSTER for example, only deals with data pertaining to the booster and 
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Figure 2: PAYLOAD INTEGRATION PROBLEM 



therefore, works independent of the payload. Also, sometimes the data 
generated by PROGRAM BOOSTER can be used in amalyses of different pay** 
loads and therefore have not to be recalculated. 

Finally, the PROGRAM approach also allows for better check-out 
and control of the data generated at several points in the process 
of a load cycle. Indeed, the user can put in his/her o%m checks if 
desired. 

PROGRAM SYSRESP respresents the hub around which the other five 
programs are centered. The purpose of PROGRAM SYSRESP is to generate 
the coupled booster/payload system response. The most important sub- 
routine called by PROGRAM SYSRESP is SUBROUTINE SRESP, which implements 
the integration scheme as outlined in Section 4 of Chapter I. The 
INPUT to PROGRAM SYSRESP consists of all the quantities necessary 
to run SUBROUTINE SRESP. The OUTPUT of PROGRAM SYSRESP is the system 
response i.e. displacements, velocities and accelerations. These 
quantities can be written on paper and tape. In particular, the pay- 
load accelerations and the booster/payload interface displacements at 
each time step are written on tape so that they can be used in PROGRAM 
LOADS for the calculation of member loads. 

Much of the INPUT to PROGRAM SYSRESP is not directly given and 
therefore must be created in advance. PROGRAM BOOSTER, PROGRAM PAYLOAD 
and PROGRAM FORCE were composed to serve this purpose. PROGRAM BOOSTER 
generates all booster data necessary to run PROGRAM SYSRESP. The 
booster organization can use this PROGRAM independently of any other 
organization. Enough subroutines were developed so that all booster 
data can be generated starting with the free mass and stiffness matrices 
[Mg] and [Kg] and the interface restrained modes and frequencies and 

[(jjg]. It is reasoncUsle to expect that these INPUT quantities are 
available. If not, the user is expected to provide this information be- 
fore running PROGRAM BOOSTER. It would not be wise to "can" the con- 

“B “2 

struction of [Mg] , [Kg] , [$^^1 and [ug] because of the multitude of ways 
these quantities can be generated. Furthermore, PROGRAM BOOSTER 
contains a number of "flags" which allow for user flexibility of INPUT. 
Indeed, often times certain quantities are already available and need 
not to be regenerated. The same is true for PROGRAM PAYLOAD whi .n is 
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very similar to PROGRAM BOOSTER except that it generates payload quan** 
titles necessary to run PROGRAM SYSRESP. In addition* it also generates 
parts of load transformations if desired. Again* much flexibility is 
possible depending on the case at hand. PROGRAM INTFACE collects sc:>e 
of the data generated by PROGRAM BOOSTER and PROGRAM PAYLOAD and pro- 
duces quantities that involve both booster and payload data. Again* 
these quantities are needed in PROGRAM SYSRESP and PROGRAM LOADS. PRO- 
GRAM INTFACE reflects the coupling between booster and payload through 
the interface. For example* it calculates the interface modes C^^j. 
PROGRAM FORCE essentially converts tne force data into the right format 
to be used in the integration program PROGRAM SYSRESP. PROGRAM <?oRCE 
also contains a number of "flags" which allows for mwre flexibility. 
Finally* as mentioned above* PROGRAM LOADS generates member loads and 
draws on PROGRAM PAYLOAD for load transformation INPUT anu on PROGRAM 
SYSRESP for payload response INPUT. 

Each of the six PROGRAMS are independent components of the soft- 
ware package. PROGRAM BOOSTER can be used by an independent booster 
organization. Similarly* PROGRAM PAYLOAD can be used by an indepen- 
dent payload organization. PROGRAM INTFACE* PROGRAM FORCE and PROGRAM 
SYSRESP can be used by an independent integration organization while 
PROGRAM LOADS can be used by any organization that is responsible for 
loads calculations. Because each of the PROGRAMS is compatible with 
the other PROGRAMS it is also possible for one organization to use 
the entire package in sequence. 

Several versions of the software package are available. All 
of these versions make use of the FORMA library. Two versions are 
available on the CDC computer (Dense-and Partition-Logic) . One ver- 
sion is available on the VAX computer (Dense) and one version is avail- 
able on the UNIVAC system (Partition-Logic) . 

This section was intended to give the reader a general idea of 
how the software package is structured. It is not intended to be a 
detailed user guide. The Final Report of this contract will contain a 
detailed user guide as well as the actual listings of all the PROGRAMS 
and associated SUBROUTINES. At this point* the software package is 
still in a "check-out" phase and is still undergoing minor changes. 
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As we shall see in subsequent sections, the software package has already 
been used on several examples and the general structure as described 
in this section is final. 
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3. 


NUMERICAL EXAMPLES 


This section we shall discuss a simple sample problem which 
is used for the initial check-out of the software package. Further- 
more , we shall briefly present the results of a realistic defense 
payload integration analysis which was conducted using the dir- 

ect integration technique. Both of these sample problems are rela- 
tively small in terms of the nvunber of modal coordinates. However, 
results produced by standard techniques are available and allow 
us to make a comparison of both techniques. The Final Report will 
contain a large scale analysis of the Shuttle Orbiter carrying 
two payloads namely: the Space Telescope and the QMS Kit. 

The first example is depicted in Figure 3 . The booster B 
consists or 18 pipe segments. The mass of each segment is equally 
divided between the end points of the segment. If we only keep trans- 
lational dofs. then the free booster has 57 dofs and the "cantilever- 
ed" booster has NB = 54 dofs. Similarly, the payload P consists of 
7 pipe segments and NP = 21. Because there are 3 rigid body modes, 
we have a determinate interface and IF * 3. The parameters for a 
booster pipe segment are: 

E=6 . 89xlO®kN/m^ , p=2 . 77xl0^kg/m^ , Aj^^^y=l . 93xl0“^m^ , 

A . «4.17xl0“^m^,L=0.762m J =8. 325xl0"^m^ 

wing o 

Similarly, for a payload pipe segment: 

E=6 . 89xl0^kN/m^ . p=2 . 77xl0^kg/m^ , A=4 . 05xl0"^m^ 

L=0.762.-n, J =1.249xl0'^m^ 
o 

Using the above dataf a finite element model was derived for botn 
the cantilevered booster and the cantilevered payload. Solving the 
eigenvalue problem yields booster frequencies ranging form 1 Kz to 
106 Hz and payload frequencies from 1 Hz to 104 H^. In this parti- 
cular example we used zero initial conditions and applied loads to 
stations 16 and 17 in the x, y, and z-directions (444822xcos (i50t) N, 
1 = 1 , 6 ) . 
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The accuracy of the adapted Newmark-Chan-Beta routine was 
checked by comparing the response results from PROGRAM SYSRESP with 
those obtained from a fourth order Runge-Kutta (Gill modification) 
routine using the same step size. The results compare very well. 

Table I shows some of the results for a step size h = 0.001 seconds. 

It should be pointed out that a phase shift in the response was ob- 
served . This is to be expected and is inherent to the numerical 
technique used. However, this phase shift does not affect the value 
of the maximum or minimum loads in an element. It may slightly af- 
fect the time at which this maximum or minimum occurs, but this is of 
little consequence. Indeed, it is fair to say that in practice every 
numerical integration technique produces a phase shift when a time 
step is used consistent with a cut-off frequency. 

Next, we compared computer cost of the present method with that 
of the conventional approach. The conventional "full-scale" approach 
first calculates the so-called "modal modes and frequencies" (i.e. 
the system modes from Equation (10)). Then, a nximerical scheme (e.g. 
Runge-Kutta) is used to determine the response from the uncoupled sys- 
tem equations. For the present exeunple, the cost of the direct inte- 
gration routine to determine the response is less by a factor of 10 
compared to the conventional approach. In the Final Report we hope 
to include a study of how this factor relates to the size of the sys- 
tem and the number of time points and interface dofs. 

We also compared the cost of the load calculations. The improv- 
ed technique decreases the cost by a factor of 12. Again, at this 
point, it is hard to tell how this factor will change when the booster/ 
payload system represents a more realistic configuration. Also, it 
may be hard to compare cost factors for large systems from a logistic 
point of view. Indeed, payload organizations currently deliver load 
transformations which are not consistent with the present improved ap- 
proach. Therefore, it is often impossible to generate the appropri- 
ate quantities required for use in the present approach, for lack of 
certain information. It should be noted that this is only a logistics 
problem and changes could easily be accommodated. 

Table II lists some of the load results obtained for this exaunple. 
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Again, note tha^ a phase- shift is present (although not noticible at 
the interface). The maximum loads however, compare extremely %iell 
(i.e. a relative error of less than 0.1%), 

In conclusion, we can state that the entire direct integration 
routine (i.e. response plus load calculations) is less expensive by 
roughly a factor of 10 compared to the conventional technique. Note 
that the analysis was performed on both the CYBER and VAX computer 
systems. The analysis of a Defense booster/payload system was also 
completed. The model consisted of a 261 degree of freedom payload model 
and a 92 degree of freedom booster model. A cut-off frequency of 50 
Hz was chosen, resulting in a coupled booster/payload model containing 
27 cantilevered booster modes and 55 cantilevered payload modes. 

There are 6 interface degrees of freedom, which makes the interface 
determinate. The entire analysis was performed on the VAX computer 
system. 

The results were very encouraging. The accuracy of the maximum 
loads was excellent (i.e. the relative error compared to the standard 
technique was < 0.5%). The cost again, was approximatley less by 
a factor of 10. 

In conclusion, we can state that the new approach shows encourag- 
ing results which warrant its application to larger booster/payload 
systems. Indeed, plans were made to apply the present technique to 
a number of Defense Systems and also to a STS-System. We intend to 
report on the results in the Final Report of this contract. 
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FINAL ROtARKS 


In this final section we wish to collect a few relati/ely unre- 
lated observations. 

First of all, there is the question of system modes and frequo’ — 
cies. Indeed, sometimes the payload integration effort is part of 
a larger program involving the solution of problems for which knowledge 
of the booster/payload system modes is required. As stated before, 
the present payload integration method does not require the calcu- 
lation of the system modes. However, these modes and frequencies 
can always be calculated from Equation (10). Usually, a one time cal- 
culation should be sufficient and no need exists to recalculate these 
system modes for every load cycle. Furthermore, the system response 
is already known and could help to reduce the cost of solving the 
system eigenvalue problem in many cases where a significant percen- 
tage of booster modes do not couple with payload modes and vice-versa. 

Next, we would like to address a problem related to the practical 
application of payload integration methods, in particular in the case 
of the Space Transportation System. As shown in Section 4 of Chapter 
I, the payload integrator must receive from the booster organization 
a set of cantilevered booster modes. In theory, this does not present 
any problem. In oractice^ however, the payload integrator often has 
to consider several alternative sets of interface points for possible 
connection of the payload. For every different interface there is a 
corresponding different set of booster modes. For reasons of logis- 
tics and cost it is virtually impossible to generate a different set 
of modes every time the interface changes. Therefore only one set of 
cantilevered booster modes is generated containing all possible inter- 
face points. As an example, let us consider a payload that has 7 in- 
terface dofs. but could be attached to the booster at 36 different 
dofs. That means the payload integrator has to select 7 appropriate 
dofs. out of the 36. The booster organization then could generate a 
set of cantilevered booster modes about the 36 interface dofs. As 
can be easily seen this creates a problem when it comes to implement- 
ing the theory of Section 4 in Chapter I. 
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This could be construed as a disadvantage of using cantilevered 
modes for the booster. Indeed, if one would use free booster modes 
in an integration technique then this problem of changing booster modes 
does not exist. However, as we shall show in the Final Report, this 
particular problem can easily be remedied. Similarly, the problem of 
multiple payloads can also be addressed in a satisfactory manner. 

C' early, the case of superfluous interface dofs. on the booster side 
and the case of multiple payloads requires a substantial change in the 
analytical write-up as well as in the computer program code. We 
shall report on these changes in the Final Report. 

The next subject deals with the cost of the response routine 
i.e. PROGRAM SYSRESP. From Section 4 of Chapter I it follows that the 
execution time of PROGRAM SYSRESP is proportional to IF, MB, MP and 
the number of time steps necessary to integrate over the desired time 
interval. The Final Report will contain a more detailed operation 
count and an approximate formula for the estimated execution time. 

At this point it should be noted that IF, NB, and NP are important 
factors in determining the cost of the integration routine. In- 
deed, it is not difficult to see that if IF is small compared to NB 
and NP the execution time goes down considerably. Similarly, the 
smaller NB and NP the less the cost. This is one clue to a possible 
short-cut method. Indeed, in many cases the coupling between the 
booster and the payload is limited. Numerically, this means that 
some rows of B and/or P will have elements equal to zero or close to 
zero, so that that particular row is decoupled from the system there- 
by effectively reducing the size of B and P without all the zero 
and "near zero" elements. This approach was also reported on in the 
Monthly Progress Report, Issue 22, Reference [12]. The question then 
is to develop a criterion defining what a "near zero" element is. 

There are many other ideas for developing short-cut methods some of 
which are outlined in Reference [ 3 ] . 

It is our intention to include a chapter on short-cut methods 
in the Final Report. 
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COMPARISON OF FOURTH-ORDER 
RUNGE-KUTTA TECHNIQUE AND 
NEWMARK-CHAN-BETA TECHNIQUE 



- 2/746 10 - 2.746 10 
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